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.«? ! ABSTRACT 

We report on a three year spectroscopic monitoring program of the Ha emis- 
sion in the massive X-ray binary LS I +65 010 = 2S 0114+650, which consists 
of a B-supergiant and a slowly rotating X-ray pulsar. We present revised orbital 
elements that yield a period of P = 11.5983 + 0.0006 d and confirm that the orbit 
has a non-zero eccentricity e = 0.18 + 0.05. The Ha emission profile is formed in 
the base of the wind of the B-supergiant primary, and we show how this spectral 
feature varies on timescales that are probably related to the rotational period of 
the B-supergiant. We also examine the X-ray fluxes from the Rossi X-ray Timing 
Explorer All-Sky Monitor instrument, and we show that the X-ray orbital light 
curve has a maximum at periastron and a minimum at the inferior conjunction 
of the B-supergiant. We also show that the wind emission strength and the high 
energy X-ray flux appear to vary in tandem on timescales of approximately a 
year. 

Subject headings: binaries: spectroscopic^ — pulsars: individual (2S 0114+650) — 
stars: early-type — stars: individual (LS I +65 010) — supergiants — stars: winds, 
outflows 



1. Introduction 

Massive X-ray binaries consist of a mas- 
sive, luminous star with a neutron star or 
black hole companion with an orbital sepa- 
ration small enough to power an accretion- 
driven X-ray source. The mass donor is of- 
ten a rapidly rotating Be star, and the neu- 
tron star companion accretes gas from the 
slowly outflowing circumstellar disk of the 
Be star (Coe 2000). A second kind of mas- 
sive X-ray binary contains an OB super- 
giant as the mass donor, and mass trans- 
fer occurs by a Roche lobe overflow stream 
or stellar wind capture (Kaper 1998). In 
both cases the mass transfer rate can be 
time variable and lead to large variations 
in the X-ray flux. The gas in the immedi- 
ate vicinity of the mass donor will generally 
be dense and hot enough to be a source of 
Ha emission, and our goal in this paper 
is to investigate how the mass loss fluctu- 
ations at the source (as observed in Ha) 
are related to the X-ray variations, particu- 
larly as measured with the All-Sky Monitor 
(ASM) instrument aboard the Rossi X-ray 
Timing Explorer (RXTE) satellite (Levine 
et al. 1996). We present similar investiga- 
tions of these co-variations in companion 
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papers on the black hole binary Cyg X-l 
= HDE 226868 (Gies et al. 2003) and the 
Be X-ray binaries LS I +61 303 (Grund- 
strom et al. 2006a) and HDE 245770 = 
A 0535+26 and X Per (Grundstrom et al. 
2006b). 

Our subject here is the massive X-ray 
binary system 2S 0114+650 (4U 0114 + 65) 
whose optical counterpart is LS I +65 010 
(V662 Cas; HIP 6081; ALS 6517). The 
characteristics of the mass donor star in 
this system were determined in a spectro- 
scopic and photometric study by Reig et 
al. (1996). They show that the star is a 
supergiant with a spectral classification of 
Bl la. They also suggest that the Ha emis- 
sion present in the spectrum forms in the 
stellar wind of the supergiant and that the 
expected wind parameters and accretion 
rate are consistent with the observed X- 
ray properties. The orbital elements were 
first determined by Crampton, Hutchings, 
& Cowley (1985) who found that the or- 
bit has a small eccentricity and a period 
of 11.6 d. The orbital period is also found 
in the X-ray flux variations (Farrell, Sood, 
& O'Neill 2006; Wen et al. 2006) along 
with two other periodic signals: a super- 
orbital period of 30.7 d (possibly related to 
the precession of the disk surrounding the 
neutron star; Farrell et al. 2006) and the 
pulsar period of 2.7 hr (Finley, Belloni, & 
Cassinelli 1992; Hall et al. 2000; Koenigs- 
berger et al. 2003; Bonning & Falanga 
2005). This is the longest spin period of 
any known X-ray pulsar, and Li & van den 
Heuvel (1999) argue that the pulsar was 
probably born with a huge magnetic field 
(as a magnetar) in order to have spun down 
sufficiently within the lifetime of the mass 
donor star. 
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Here we outline the results of an Ha 
monitoring program on LS I +65 010 that 
spanned the years 1998 - 2000. We de- 
scribe the spectroscopy in §2 and present 
revised orbital elements in §3. We discuss 
the properties and temporal variations of 
the Ha emission in §4, and in §5 compare 
them with the X-ray variations observed 
by RXTE/ ASM. 

2. Observations 

We obtained 100 optical spectra of 
LS I +65 010 with the Kitt Peak National 
Observatory 0.9 m coude feed telescope 
during six observing runs between 1998 
August and 2000 December. The details 
of the spectrograph arrangement for each 
run are listed in Table 1. All the spec- 
tra cover the region surrounding Ha, but 
they were made using two different spec- 
tral dispersions. During the first, fifth, 
and sixth runs, we used grating B (632 
grooves mm -1 with a blaze wavelength of 
6000 A in second order) and obtained a 
resolving power of R — X/8X ~ 9500. The 
spectra from the other three runs were 
made with the lower dispersion RC181 
grating (316 grooves mm -1 with a blaze 
wavelength of 7500 A; made in first order 
with a GG495 filter to block higher or- 
ders), which yielded an average resolving 
power of R = 4100 — 4400. Exposure times 
were usually 30 minutes, and the spectra 
generally have aS/N» 100 pixel -1 in the 
continuum. 

The spectra were extracted and cali- 
brated using standard routines in IRAF 5 , 

6 IRAF is distributed by the National Optical As- 
tronomy Observatory which is operated by the 
Association of Universities for Research in Astron- 



and then each continuum rectified spec- 
trum was transformed onto a uniform 
heliocentric wavelength grid for analysis. 
The removal of atmospheric lines was done 
by creating a library of spectra from each 
run of the rapidly rotating A-star ( Aql, 
removing the broad stellar features from 
these, and then dividing each target spec- 
trum by the modified atmospheric spec- 
trum that most closely matched the target 
spectrum in a selected region dominated 
by atmospheric absorptions. 

3. Orbital Radial Velocity Curve 

The spectra from different runs cover 
different portions of the red spectrum (Ta- 
ble 1), but all record both the Ha emis- 
sion line and the He I A6678 absorption 
line. Thus, we will concentrate on these 
common features in the analysis. We mea- 
sured radial velocities by cross-correlating 
each spectrum with a standard template 
that was formed by averaging all the spec- 
tra from the fifth run (2000 October). The 
spectral range for the cross-correlation was 
generally limited to the immediate vicinity 
of the He I A6678 line, but we also included 
the region surrounding He I A7065 for the 
spectra from the final two runs (2000 Octo- 
ber and December). The positions of these 
two features in the template spectrum were 
measured by parabolic fits to the lower 
portions of the profiles, and we added their 
mean radial velocity, —52.8 ± 0.4 km s -1 , 
to each relative velocity to place them on 
an absolute scale. We note that there are 
known line-to-line velocity differences in 
the spectrum that are related to line for- 

omy Inc., under cooperative agreement with the 
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mation at different depths in the expand- 
ing atmosphere (Crampton et al. 1985; 
Koenigsberger et al. 2003), and velocities 
based on the strong He I lines are probably 
systematically more negative than those 
from lines formed deeper in the atmo- 
sphere. We also made cross-correlation 
measurements of the strong interstellar line 
at 6613 A that we used to monitor the ra- 
dial velocity stability of our spectra. Our 
results are summarized in Table 2 (which 
appears in full in the electronic version) 
that lists the mid-point heliocentric Ju- 
lian date of observation, the orbital phase, 
radial velocity, observed minus calculated 
residual from the fit, the measured veloc- 
ity difference of the 6613 A interstellar line 
from its position in the template spectrum, 
and the equivalent width of the Ha emis- 
sion (§4). The typical stellar radial veloc- 
ity errors are about ±1.7 km s _1 for the 
high dispersion spectra and about twice 
that value for the lower dispersion spec- 
tra (based upon the differences of closely 
spaced pairs of observations). 

We determined orbital elements from 
these velocities using the non-linear least- 
squares fitting program of Morbey & 
Brosterhus (1974). We initially deter- 
mined elements for both the new data and 
the archival measurements from Crampton 
et al. (1985), and while most of the ele- 
ments were in good agreement, we found 
the systemic velocity 7 for the new set was 
2.5 km s _1 lower than that based on the 
data from Crampton et al. (1985). This is 
not surprising given the existence of sys- 
tematic line-to-line velocity differences and 
that Crampton et al. (1985) measured a 
different set of lines in the blue. Conse- 
quently, we subtracted this value from all 



the archival velocities, and then made a 
new fit of the merged velocities. The re- 
vised elements are presented in Table 3, 
and the radial velocity curve is shown in 
Figure 1. The revised orbital elements are 
in reasonable agreement with the original 
values from Crampton et al. (1985) that 
are also listed in Table 3. The main differ- 
ences are an order of magnitude improve- 
ment in the accuracy of the orbital period 
and a modest decrease in the orbital semi- 
amplitude K. 

Crampton et al. (1985) were unable to 
decide between circular and elliptical solu- 
tions, but it is now clear that a non-zero 
eccentricity is required. We made fits as- 
suming both circular and elliptical param- 
eters, and the improvement in the resid- 
uals in the elliptical case would only oc- 
cur by random chance with a probability 
of 0.8% according to the test of Lucy & 
Sweeney (1971). Thus, the eccentricity is 
significantly different from zero, and we 
report only the elliptical solution in Ta- 
ble 3. The residuals from the fit (rms = 
6.4 km s _1 ) are larger than those associ- 
ated with the measurement errors (1.7 - 
3.4 km s -1 ), and, in fact, we zero-weighted 
three measurements that had exceptionally 
large deviations (Fig. 1). We suspect that 
this scatter is due to the microvariablility 
inherent to the atmospheres of early-type 
supergiants (Gies & Bolton 1986). 

4. Ha Variations 

The Ha feature has always appeared 
in emission in past observations (Reig et 
al. 1996; Liu & Hang 1999; Koenigsberger 
et al. 2003) and emission was consistently 
present in our spectra as well. However, 
the shape and strength of the profile var- 



4 



ied on a timescale of days. We show in 
Figure 2 two example time sequences from 
our best sampled runs. A red-shifted emis- 
sion peak is generally present and in about 
half of the spectra a blue-shifted absorp- 
tion trough is also visible. There are indi- 
cations in several cases of absorption fea- 
tures that appear to migrate bluewards 
to velocities of ps —400 km s _1 . These 
kinds of variations are common among B- 
supergiants (Fullerton et al. 1997; Markova 
& Valchev 2000) and probably result from 
transient structures in the wind. Simi- 
lar kinds of blueward absorption variations 
were also observed in the He I A5876 line in 
the spectrum of LS I +65 010 by Koenigs- 
berger et al. (2003). Liu & Hang (1999) 
interpreted episodes of double-peaked Ha 
emission as evidence of a second emission 
component from the vicinity of the neutron 
star companion, but our time sequences 
suggest that they are instead the result of 
varying stellar wind absorption. 

We measured emission equivalent widths 
for Ha by making a simple numerical inte- 
gration of line flux above the continuum, 
and our results are listed in the final col- 
umn of Table 2. The typical measurement 
errors are approximately ±0.05 A for the 
higher dispersion spectra and about twice 
that for the lower dispersion spectra. The 
resulting equivalents widths fall within the 
range reported in the past (Reig et al. 1996; 
Liu & Hang 1999). We searched for any pe- 
riodic signals in the time variations of the 
Ha equivalent widths using the discrete 
Fourier transform and CLEAN algorithm 7 
(Roberts, Lehar, & Dreher 1987). These 
results are summarized for each observ- 



7 Writtcn in the Interactive Data Language by A. 
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ing season and for the entire data set in 
Figure 3. Most of the temporal variabil- 
ity appears to occur on timescales around 
10 d (frequency of 0.1 cycles d _1 ). We 
found no compelling evidence of periodic 
variability at either the orbital (11.6 d) 
or super-orbital period (30.7 d), and the 
la upper limits for such cyclic variations 
are < 4% and < 10% for these two pe- 
riods, respectively (based upon sinusoidal 
fits at these fixed periods). The strongest 
signal in the power spectrum (among a 
number of candidate periods in the range 
between 6 and 12 d) occurs at a period of 
9.72 ± 0.02 d. The false alarm probabil- 
ity that this peak is due to noise is small, 
p iVexp(— zq) = 0.0002 where z = 15.6 
is the ratio of the power peak to the vari- 
ance and iV = 1000 is the number of quasi- 
independent frequencies tested over the 
range in physically interesting timescales 
between and 1 cycle d _1 (Scargle 1982). 
Although Ha variations occurred during 
each run on such timescales, there were 
large differences in the mean equivalent 
width between runs and in the detailed 
profile shape variations within runs that 
indicate the importance of other timescales 
as well. We suggest that this 10 d timescale 
represents a characteristic recurrence time 
for ephemeral structures to appear in a 
wind that is subject to longer term vari- 
ations in mass loss rate. If the structures 
are azimuthally distributed with respect to 
the star's spin axis, then this characteristic 
time may correspond to the stellar rotation 
period (see the case of the B-supergiant 
HD 64760 discussed by Fullerton et al. 
1997 and the models for rotational mod- 
ulation presented by Harries 2000). We 
can estimate the rotational period very 
approximately using the projected rota- 
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tional velocity Vsmi = 96 ± 20 km s^ 1 
and radius estimate R/R & = 37 + 15 (Reig 
et al. 1996) plus the probable inclination 
range (see Fig. 6 below). The timescale for 
the Ha variations falls within the derived 
range of rotational period of 9 to 15 d. 

5. X-ray Variations 

The 2S 0114+650 system is one of some 
hundred targets that have been under con- 
tinuous surveillance since 1996 March by 
the All-Sky Monitor on the Rossi X-ray 
Timing Explorer (Levine et al. 1996). This 
instrument records the X-ray flux across 
three energy bands in 90 s observation seg- 
ments. We extracted the flux measure- 
ments from 1996 March through 2006 May 
from the quick-look results provided by 
the RXTE/ ASM team 8 . The source is 
relatively weak and there is a significant 
amount of scatter among the individual 
measurements. However, since there are 
currently over 60,000 measurements avail- 
able for 2S 0114+650, we can achieve a sig- 
nificant increase in S/N by averaging the 
fluxes into temporal or orbital-phase bins 
of interest. 

We first consider the long term varia- 
tions in X-ray flux. We show in the up- 
per panel of Figure 4 the high energy band 
fluxes for the duration of the RXTE mis- 
sion binned into intervals equivalent to 12 
binary orbits (139 d). Each symbol shows 
the mean and standard deviation of the 
mean of the flux in the time bin. These 
means and errors were taken from mea- 
surements made only with the Scanning 
Shadow Camera (SSC) 2 and SSC3, since 
there is some indication of a long term 

8 http://xte. mit.edu/ 



trend of uncertain origin in the results from 
SSC1 for this particular source (possibly 
related to the calibration of gain changes 
in the proportional counter associated with 
SSC1). The lower panel illustrates the Ha 
equivalent width measurements over the 
same dates (but with much more restricted 
time coverage). There is a significant range 
in observed equivalent width during each 
run (§4), but it appears that mean strength 
was largest in 1998 (~ HJD 2,451,060), 
smallest in 1999 (~ HJD 2,451,460), and 
intermediate in 2000 (~ HJD 2,451,860). 
The binned X-ray fluxes in the top panel 
also show a similar kind of variation, which 
is suggestive that the time-averaged Ha 
emission and X-fluxes appear to vary to- 
gether. Taken at face value, this indi- 
cates that long term variations in the B- 
supergiant mass loss rate are reflected in 
the X-ray accretion flux. A similar result 
was found for the massive X-ray binary and 
microquasar LS 5039 (Reig et al. 2003; Mc- 
Swain et al. 2004). 

Next we binned the ASM fluxes accord- 
ing to orbital phase from our spectroscopic 
solution (Table 3), and the light curves for 
the three energy bands are plotted in Fig- 
ure 5. The orbital light curve is best de- 
fined in the 5-12 keV band where the 
observed counts are highest (top panel of 
Fig. 5; compare with the earlier light curve 
presented in Fig. 1 of Hall et al. 2000). We 
see that there is a clear maximum when 
the stars are closest near periastron (phase 
0.0). At that point in the orbit, the neu- 
tron star will encounter the densest and 
slowest portions of the supergiant's wind, 
and the wind accretion rate will attain 
a maximum (Lamers, van den Heuvel, & 
Petterson 1976). The highest 5-12 keV 
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flux average actually occurs for the bin cen- 
tered at phase 0.08, corresponding to a 
time some 22 hours after periastron. This 
may indicate the interval between accre- 
tion and heating to X-ray emitting tem- 
peratures or it may represent the transit 
time for a wind density enhancement cre- 
ated by tidal effects at periastron to reach 
the vicinity of the neutron star. 

The minimum in the X-ray light curve 
appears near phase 0.63, which is very 
close to the predicted phase of supergiant 
inferior conjunction at phase 0.65 ± 0.06. 
The X-ray flux does not vanish at this 
minimum, so we suspect that we are ob- 
serving an atmospheric rather than a to- 
tal eclipse (Hall et al. 2000). For example, 
Wen et al. (1999) show that in the case 
of Cyg X-l, the decrease in X-ray flux ob- 
served when the supergiant is in the fore- 
ground is well explained by the attenuation 
of the X-ray flux by the wind of the super- 
giant. We suspect that the same explana- 
tion applies to the case of the LS I +65 010 
= 2S 0114+650 system since Hall et al. 
(2000) find evidence of an increase in col- 
umn density near X-ray minimum. 

The kind of analysis made by Wen et al. 
(1999) of the wind attenuation func- 
tion of orbital phase would also be prof- 
itable in this case and would presumably 
help place limits on the orbital inclina- 
tion. We show in Figure 6 the constraints 
on the masses of the two components in 
the mass plane diagram. The solid lines 
show the mass relations for given values 
of orbital inclination that are derived from 
the mass function (/(m) given in Table 3). 
Reig et al. (1996) estimate that the su- 
pergiant mass and radius are 16 + 5 M 
and 37+15 R & , respectively, based upon 



their analysis of the optical and near-IR 
spectrum. Note that this supergiant ra- 
dius is much smaller than the star's Roche 
radius for the full range of acceptable in- 
clination and mass ratio. The limiting line 
for the absence of full eclipses is shown by 
the dashed line in Figure 6 (derived from 
the radius given above and the projected 
semimajor axis ctisini listed in Table 3). 
Finally, we can place a weak constraint on 
a minimum inclination angle of i = 20° 
by assuming that the supergiant rotates no 
faster than the critical rate (based on the 
given radius and observed projected rota- 
tional velocity V sini = 96 ± 20 km s -1 ; 
Reig et al. 1996). This is a weak con- 
straint because it is based upon the as- 
sumption that the supergiant 's spin incli- 
nation is the same as the orbital incli- 
nation, but these inclinations could dif- 
fer if the neutron star suffered an asym- 
metric kick at the time of the supernova 
explosion (Brandt & Podsiadlowski 1995). 
The shaded region in Figure 6 shows the 
probable B-supergiant range together with 
the observed range in neutron star mass 
(van Kerkwijk, van Paradijs, & Zuiderwijk 
1995), and a plus sign marks the nominal 
position for a Chandrasekhar mass neutron 
star (1.4 M ). This best estimate solution 
has an inclination of i = 45°, just below 
eclipse limit near i = 50°, and this orien- 
tation is probably consistent with the at- 
mospheric eclipse seen in the X-ray light 
curve. A realistic model of the X-ray at- 
tenuation by the wind of the B-supergiant 
would provide a reliable estimate of the 
inclination and better constraints on the 
masses of both stars. 
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Table 1 

Journal of Optical Spectroscopy 



Dates 


Spectral Range 


Resolving Power 


Number 


(HJD-2,400,000) 


(A) 


(A/AA) 


of Spectra 


51053 - 51065 


6313 - 6978 


9530 


23 


51421 - 51429 


5397 - 6735 


4050 


11 


51464 - 51469 


5400 - 6736 


4100 


14 


51491 - 51497 


5545 - 6881 


4400 


14 


51817-51830 


6440 - 7105 


9500 


16 


51890 - 51901 


6443 - 7108 


9500 


22 
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Table 2 

Radial Velocity and Equivalent Width Measurements 



Date 


Orbital 


V r 


(O - o 


AVism 


W\(H.a) 


(HJD-2,400,000) 


Phase 


(km s _1 ) 


(km s — ) 


(km s~ ) 


(A) 


51053.803 


0.481 


-70.8 


-4.5 


-1.8 


-2.05 


51053.824 


0.483 


-66.8 


-0.5 


-0.4 


-2.01 


51055.810 


0.655 


-55.1 


1.1 


0.6 


-1.98 


51055.832 


0.656 


-56.3 


-0.2 


0.7 


-2.10 


51055.917 


0.664 


-53.9 


1.8 


3.9 


-2.02 


51056.826 


0.742 


-56.1 


-6.0 


2.4 


-1.70 


51056.850 


0.744 


-54.4 


-4.5 


2.3 


-1.68 


51056.871 


0.746 


-57.0 


-7.3 


2.2 


-1.74 


51057.810 


0.827 


-50.5 


-5.8 


1.0 


-1.53 


51057.831 


0.829 


-49.7 


-5.2 


2.0 


-1.59 


51057.852 


0.831 


-49.7 


-5.3 


2.1 


-1.66 


51058.801 


0.912 


-47.0 


-4.4 


3.6 


-1.15 


51058.822 


0.914 


-49.8 


-7.1 


2.8 


-1.16 


51058.844 


0.916 


-49.2 


-6.5 


2.4 


-1.28 


51061.800 


0.171 


-61.0 


5.0 


-2.0 


-1.77 


51061.821 


0.173 


-57.6 


8.6 


-0.1 


-1.77 


51061.843 


0.175 


-59.2 


7.1 


2.3 


-1.73 


51063.887 


0.351 


-73.2 


-2.9 


2.0 


-1.11 


51063.908 


0.353 


-71.4 


-1.0 


1.7 


-1.13 


51063.930 


0.355 


-74.6 


-4.3 


0.5 


-1.17 


51065.811 


0.517 


-71.5 


-6.9 


0.5 


-1.54 


51065.835 


0.519 


-66.0 


-1.5 


-0.5 


-1.63 


51065.856 


0.521 


-64.4 


0.0 


0.7 


-1.75 


51421.915 


0.220 


-69.0 


-0.2 


-3.2 


-0.98 


51423.896 


0.391 


-75.5 


-5.9 


-3.2 


-1.30 


51425.898 


0.563 


-68.8 


-6.7 


-6.6 


-1.27 


51425.919 


0.565 


-67.4 


-5.4 


-7.9 


-1.28 


51426.879 


0.648 


-61.4 


-4.7 


-0.4 


-0.82 


51426.901 


0.650 


-60.9 


-4.3 


-0.2 


-0.85 


51427.879 


0.734 


-56.3 


-5.7 


1.2 


-0.86 


51428.838 


0.817 


-48.2 


-3.1 


-1.8 


-0.79 


51428.860 


0.819 


-48.7 


-3.6 


-1.1 


-0.82 


51429.836 


0.903 


-46.1 


-3.5 


-1.5 


-0.93 


51429.857 


0.905 


-46.4 


-3.8 


1.0 


-0.90 


51464.760 


0.914 


-44.5 


-1.8 


-0.3 


-1.43 


51464.783 


0.916 


-41.9 


0.9 


1.0 


-1.17 


51465.814 


0.005 


-57.9 


-9.7 


-3.8 


-1.57 


51465.835 


0.007 


-56.0 


-7.6 


-5.6 


-1.63 


51465.858 


0.009 


-57.5 


-8.8 


-4.0 


-1.50 


51466.772 


0.087 


-56.0 


1.9 


-8.4 


-1.28 


51466.794 


0.089 


-57.1 


1.1 


-5.2 


-1.26 


51466.816 


0.091 


-59.0 


-0.6 


-3.4 


-1.35 


51467.855 


0.181 


-72.7 


-5.9 


-0.4 


-1.03 


51467.876 


0.183 


-70.7 


-3.9 


-2.2 


-1.19 


51468.812 


0.263 


-70.9 


-0.7 


1.3 


-1.46 


51468.834 


0.265 


-65.4 


4.8 


0.8 


-1.49 


51469.829 


0.351 


-67.6 


2.8 


-6.1 


-1.28 


51469.850 


0.353 


-70.5 


-0.2 


-2.1 


-0.72 
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Table 2 — Continued 



Date 


Orbital 


V r 


(O-C) 




W x (Ua) 


(HJD-2,400,000) 


Phase 


(km s- 1 ) 


(km s- 1 ) 


(km s- 1 ) 


(A) 


51491.753 


0.241 


-56.9 


12.7 


-2.9 


-1.56 


51491.775 


0.243 


-58.0 


11.7 




-1.64 


51492.713 


0.324 


-64.2 


6.5 


3.0 


-1.70 


51492.734 


0.326 


-69.8 


0.8 


-0.8 


-1.66 


51493.696 


0.409 


-71.9 


-2.8 


3.0 


-1.46 


51493.717 


0.411 


-73.0 


-3.9 


3.4 


-1.43 


51494.708 


0.496 


-60.3 


5.4 


-2.3 


-1.14 


51494.730 


0.498 


-58.9 


6.7 


-4.5 


-1.22 


51495.765 


0.587 


-62.5 


-1.9 


1.1 


-0.91 


51495.787 


0.589 


-58.9 


1.6 


1.1 


-0.89 


51496.757 


0.673 


-66.3 


-11.3 


-4.3 


-0.95 


51496.778 


0.675 


-66.2 


-11.3 


-4.1 


-0.46 


51497.718 


0.756 


-55.6 


-6.5 


-2.5 


-0.87 


51497.739 


0.757 


-56.4 


-7.4 


-1.6 


-0.81 


51817.702 


0.344 


-69.6 


0.8 


-1.0 


-1.13 


51817.739 


0.348 


-70.3 


0.2 


-1.9 


-1.20 


51818.722 


0.432 


-64.1 


4.2 


-1.0 


-1.16 


51818.746 


0.434 


-68.3 


-0.0 


0.3 


-1.14 


51819.693 


0.516 


-59.7 


5.0 


1.0 


-1.10 


51819.716 


0.518 


-59.9 


4.7 


-2.4 


-1.02 


51820.735 


0.606 


-56.3 


3.2 


-0.0 


-1.21 


51821.698 


0.689 


-52.1 


1.7 


-0.6 


-1.20 


51821.720 


0.691 


-52.5 


1.2 


-1.5 


-1.24 


51822.710 


0.776 


-24.2 a 


23.5 


1.3 


-1.21 


51823.692 


0.861 


-22.2 a 


21.1 


0.2 


-1.65 


51823.713 


0.863 


-21. a 


22.2 


0.7 


-1.56 


51824.687 


0.947 


-42.9 


0.8 


3.7 


-1.32 


51824.708 


0.949 


-43.8 


0.1 


4.1 


-1.30 


51830.713 


0.466 


-64.5 


2.5 


-1.2 


-1.89 


51830.735 


0.468 


-65.0 


1.9 


0.2 


-1.75 


51890.622 


0.632 


-65.4 


-7.6 


-2.0 


-1.81 


51890.643 


0.633 


-66.2 


-8.5 


0.8 


-1.86 


51892.613 


0.803 


-45.3 


0.6 


-2.9 


-1.86 


51892.635 


0.805 


-45.1 


0.8 


-1.4 


-1.76 


51893.642 


0.892 


-43.0 


-0.4 


-0.3 


-1.39 


51893.664 


0.894 


-43.0 


-0.4 


0.0 


-1.34 


51894.640 


0.978 


-51.1 


-5.4 


-2.0 


-1.42 


51894.661 


0.980 


-51.5 


-5.5 


-0.5 


-1.33 


51895.620 


0.063 


-60.5 


-5.5 


-0.7 


-1.51 


51895.642 


0.064 


-59.4 


-4.2 


-0.1 


-1.49 


51896.662 


0.152 


-70.3 


-5.7 


-0.5 


-1.07 


51896.683 


0.154 


-70.9 


-6.2 


0.5 


-1.24 


51897.664 


0.239 


-67.1 


2.5 


-1.5 


-1.15 


51897.685 


0.241 


-69.6 


0.0 


-0.6 


-1.10 


51898.667 


0.325 


-72.0 


-1.4 


-1.6 


-1.47 


51898.688 


0.327 


-72.6 


-2.0 


-1.4 


-1.50 


51899.670 


0.412 


-73.0 


-4.0 


0.3 


-1.70 


51899.691 


0.413 


-74.5 


-5.5 


-1.0 


-1.74 
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Table 2 — Continued 



Date Orbital V r (O - C) AVism W x (Ua) 

(HJD-2,400,000) Phase (km s" 1 ) (km s" 1 ) (km s" 1 ) (A) 



51900.664 0.497 -55.9 9.7 1.0 -1.67 

51900.685 0.499 -54.2 11.3 1.6 -1.64 

51901.649 0.582 -44.9 16.0 1.2 -1.31 

51901.671 0.584 -48.4 12.5 -1.1 -1.45 



a Assigned zero weight in the orbital solution. 
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Table 3 
Orbital Elements 



Element Crampton et al. (1985) This Work 



P (days) 11.588 ±0.003 11.5983 ± 0.0006 

T (HJD-2,400,000) 44134.9 ± 0.7 51825.3 ± 0.5 

e 0.16 ±0.07 0.18 ±0.05 

cu (deg) 11 ± 11 51 ± 17 

K (km s" 1 ) 17 ± 1 14.0 ±0.7 

7 (km s" 1 ) -57 ±1 -58.2 ± 0.5 

/(to) (M ) 0.006 ± 0.001 0.0032 ± 0.0005 

a lS mi(R Q ) 3.9 ± 0.2 3.16 ±0.16 

r.m.s. (km s" 1 ) . . . 6.8 6.4 
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Fig. 1. — Calculated radial velocity curve 
(solid line) for LS I +65 010. The mea- 
sured radial velocities are shown as filled 
circles. Open circles mark the measure- 
ments that were assigned zero weight (§3), 
and plus signs indicate the measurements 
from Crampton et al. (1985). 
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Fig. 2. — Temporal variations in the Ha 
profile. The left panel shows the pro- 
files from 1998 August while the right 
panel shows the variations in 2000 De- 
cember. The vertical bar at upper right 
indicates the spectral intensity scale rel- 
ative to a unit continuum. The spectra 
were smoothed with a Gaussian filter with 
FWHM = 30 km s _1 . 
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Fig. 3. — Periodograms of the Ha equiv- 
alent width variations for data from 1998, 
1999, 2000, and all years (top to bottom, 
with offsets of 0.6, 0.45, 0.3, and 0.15, 
respectively). The lower plot shows the 
CLEANed periodogram (offset by 0.05) de- 
rived from a deconvolution of the full sam- 
ple periodogram. 
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Fig. 4. — The long term variations in 
the RXTE/ ASM high energy flux and Ha 
equivalent width. The top panel shows the 
ASM counts binned into intervals of 139 d 
(12 orbital periods) while the lower panel 
illustrates the Ha emission strength. The 
plus signs indicate our measurements while 
the single x sign is based on the observa- 
tion from Koenigsberger et al. (2003). 
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Fig. 5.— The RXTE/ ASM fluxes binned 
according to orbital phase (phase 0.0 for 
periastron) . The panels show the bin mean 
and standard deviation of the mean of 
ASM counts for each of the three energy 
bands. 
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5F 




M B (SOLAR MASSES) 



Fig. 6. — The relations between the super- 
giant mass M B and neutron star mass M x 
for LS I +65 010. Lines of constant or- 
bital inclination are constructed from the 
mass function and are labeled at the right 
hand side. The dashed line indicates the 
upper limit on inclination and the lower 
limit on X-ray source mass based upon the 
lack of eclipses. The shaded region shows 
the probable solution space with the plus 
sign marking the favored solution. 



19 



